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^ I New properties of the muhiplicity distributions predicted by higher order QCD and their 

physical origin are discussed briefly. Several studies which can be performed at HERA 
are proposed. 
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1 Definitions 



Multiplicity distributions (MD) of particles produced in high energy collisions are the most 
typical and widely discussed characteristic of the interaction dynamics. In a condensed 
form MD provide information about the fluctuations of energy spent for multiple particle 
production during a collision. 

The goal of the present paper is to review briefly the new features of the multiplicity 
distributions predicted by higher order QCD. 

There are two complementary ways of dealing with multiplicity fluctuations : 

- studying the distribution P„ = Onjo which is the number of produced particles per 
event, or 

- measuring the inclusive multiplicity correlators. 

In practice, one uses often the normalized factorial moments Fq and cumulants Kq 
(for review see 0) defined as 

F, = t - - « + = ^-^■■^-"^^^ , (1) 

n=0 \"/ 

q-1 

Kq = Fq~Y. C^-lKq-mFm. (2) 

m=l 

Here C!" = —n-^ — rr are the binomial coefficients and Fn = Fi = Ki =1. 

<? m\(q—my. u i ± 

These moments have an important advantage over the original moments 0. The 
average shown in (1) implies mean value of the corresponding expressions over the available 
set of experimental events. In experiment this averaging takes into account both statistical 
and dynamical effects. If one assumes that random fluctuations due to limited number of 
detected particles are described by the Poissonian distribution, then the total average of 
the factorial moments is equivalent to the dynamical average of usual moments . 

In the Feynman diagram's language, Fg corresponds to the set of all graphs while the 
cumulants Kg describe the connected graphs only. The cumulants provide the knowledge 
about the "true" correlations , non-reducible to the product of the correlations of lower 
orders. At asymptotic energies the normalized factorial moments (as well as the ordinary 
ones) do not depend on energy and are the functions of their rank only. The higher the 
rank of the moment is the more sensitive Fg and Kg are to the "tail" of MD at large n. 
The steeper decrease of the distribution at large n leads to smaller values of the high rank 
factorial moments . 

In a theoretical analysis instead of studying of the numerical series P„ it is more 
convenient to analyse the function "generating" it, namely the generating function (GF). 
Fg and Kg are easily calculated if the generating function G{u) is known 
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Thus, the knowledge of GF gives us a possibihty to calculate both the multiplicity dis- 
tribution and cumulant and factorial moments i.e. (3)- (6) demonstrate mathematical 
equivalence the description of MD by functions P„, Pg and Kg . In |^ it has been pro- 
posed to use the ratio of cumulant to factorial moments Hg = Kg/ Fg which behaves 
in a qualitatively different way for various distributions and is more sensitive to specific 
features of P„ which are invisible when just plotted P„ or even Fg (see Sec. 3). 



2 Some properties of the multiplicity distributions 

In pre-QCD time Koba, Nielsen and Olesen published the paper |^ with a hypothesis 
about the scaling properties of the multiplicity distributions at asymptotic energies (the 
KNO scaling). If z is the scaled multiplicity z = n/{n), then the KNO scaling implies a 
universal form 

^Piz) = (n) P„ 

for the multiplicity distribution. During last 30 years the KNO-like behavior of MD was 
experimentally confirmed in various types of high energy particle production processes 
except the data on proton-antiproton interactions at the highest energies y/s = 546 and 
900 GeV obtained by UA5 collaboration |[ in CERN. 
The negative binomial distribution (NBD) 



G{u) = (l - 



u{n) 



n\{k-l)\ \l + {n)/kj V k 
_ (A: + l)...(fc + g-l) _ (g-1)! 

' (fc + l)...(fc + g-l) 

is another example of the distribution which is in good agreement with experimental 
data in full phase space and in smaller phase space domains. NBD depends on two 
parameters, the average multiplicity (n) and a positive parameter k describing the shape 
of the distribution. Here we will mention only two classes of mechanisms proposed to 
generate NBD, (partial) stimulated emission || and cascading 

One feature of Hg for NBD is that it always positive and tends to zero with a q~'^ 
behaviour at high ranks. For the Poisson distribution Hg is identically equal to zero 
(except for Hi = l) 



3 



3 What does QCD tell us about the multiplicity dis- 
tributions ? 



The KNO hypothesis was strongly supported by QCD when the equations for generat- 
ing function were solved in the so-called double logarithmic approximation (DLA). DLA 
happens to be too crude however for making reasonable predictions even for asymptotic 
energies: the predicted KNO shape of the distribution appeared to be much wider than 
the experimental one. On the qualitative level, DLA can be thought to overestimate cas- 
cading processes, ignoring completely energy-momentum conservation since the energy of 
the radiating particles remains unchanged after a soft gluon emission. Therefore DLA 
apparently overestimates the gluon multiplicity because the parton characteristic energy 
is higher and the parton multiplicate more actively. Taking into account higher order per- 
turbative corrections leads to a more accurate control over the parton splitting processes 
and energy conservation. 

Such an approach has been realized (see 0, [|T0[) in the framework of the modified lead- 
ing logarithmic approximation (MLLA) by a generalization of the standard LLA scheme 
following the logic of the famous Dokshitzer-Gribov-Lipatov-Altarelli-Parisi approach and 
including the exact angular ordering (AO) (instead of the strong AO within DLA). Thus 
the system of the MLLA integro-differential equations for the quark and gluon GF has 
been derived. 

A recent series of publications [||, ||rT|-||14| was devoted to solving of these equation 
in the case of e"'"e~-collisions with account of different next-to-next-to leading (NNL) ef- 
fects. Corresponding corrections can be looked upon as being due to a more accurate 
account of energy conservation in the course of parton splitting. For example, the approxi- 
mation used in [0| allowed in the framework of gluodynamics the derivation of analytical 
expressions for the asymptotic behaviour of factorial moments and the KNO function, 
which are in better agreement with the data, by reducing substantially the width of the 
theoretical distribution. Cumulant and factorial moments of the multiplicity distributions 
in the perturbative gluodynamics have been calculated in [Q, [|l^. Accounting for the 
degrees of freedom associated with quarks |jl3[ does not change the essential qualitative 
features of Fq , Kq and influence only weakly Hq . The exact solutions of the QCD 
equations for quark and gluon GF are obtained for the case of fixed coupling in [|1^ . 

The ratio Hq is more sensitive to the form of P„ at large n than Fq (see Fig. 1). It 
was shown in |TB| that the predictions of Fq shown in Fig. la, have qualitatively the same 
behaviour and are very close to each other for q < 10. However Hq (Fig. lb) demonstrate 
much stronger sensitivity to the assumptions used. The most typical feature of the ratio 
Hq predicted by QCD [|12| is its quasi-oscillating form with a changing sign (Fig. 2). Such 
an oscillating behavior of Hq is a specific property of higher order QCD. Less complete 
account of nonlinearities in the equation for GF leads |Q, ||13[ only to one minimum with 
a very small value of Hq (the solid line in Fig. lb). 

The results of |Q have initiated a search for the peculiarities of Hq from the ex- 
perimental data. According to the Hq measurements from multiplicity distributions in 
e"'"e~-annihilation in the energy range from 22 to 91 GeV, and in /i/i-coUisions, in the 
energy range from 24 to 900 GeV, made in [0, its behaviour corresponds to the predic- 
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Fig.l. a) The factorial moments for dif- 
ferent QCD distributions ([^ - solid line; 
Tl] - dotted line) and for NBD with 
/c= 7. 6 (dashed line), b) The ratio Hq for 
the same distributions as in a). 
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Fig. 2. The ratio Hq predicted by QCD 
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tions of higher order QCD. A few examples are presented in Fig. 3. It is a surprise for 
us that the theoretical results obtained for hard processes at asymptotic energies, are in 
a qualitative agreement with experimental data at low and high energies both for e~^e~ 
processes and soft hadronic collisions. 

The behaviour of Hg for NBD shown in Fig 1, where Hq falls monotone but always 
positive, tending to zero at large ranks q, is not compatible to results shown in Fig 3. 
Therefore, despite the fact that NBD fits experimental MD very well it is not appropriate 
for the complete of description MD in particle production processes as claimed in and 
| T6|] . However, as will be seen from Sec. 5, after modifications NBD is able to generate the 
oscillating Hg as well. 



4 Monte Carlo Generators 

All Monte Carlo (MC) generators for high energy physics and, in particular, those 
which simulate deep inelastic scattering (DIS) ||18| are based on the leading logarihm 
(LL) picture with two body parton splitting a ^ d + c. However, as one mentioned in 
the previous section, higher orders in the perturbative QCD are necessary for a proper 
description of mult ipro duct ion at high energies. 

At present this can only be achieved in the generators through approximate methods 
implemented in different QCD cascades e.g. the Lund parton shower (PS) the color 
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Fig. 3. Experimental data fl^ on Hg for 
a)-d) e+e" (7^=29, 34.8, 43.8, 91 GeV) 
and e)-h) hh (v^=62.2, 200, 546, 900 
GeV) collisions. Lines are to guide the 
eye. 
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Fig. 4. The ratio Hg due to the QCD 
MC codes: a) JETSET 7.3, e+e", y/s= 
91 GeV; b) ARIADNE 4.4, e+e', ^= 91 
GeV; c) PYTHIA 5.5, e'p, ^s= 314 GeV. 
Lines are to guide the eye. 



dipole model (CDM)|^. The LLA used in PS and CDM does not give a proper treatment 
of hard emissions. A method was developed to let a single hard emission to be controlled 
by the exact 0{as) or 0{a1) QCD matrix elements and then modelling subsequent radi- 
ation using the PS technique. 

One can ask a question: are the above-mentioned improvements of the MC models 
enough for a proper description of Hq ? The answer seems to us obvious: since LLA is the 
base of PS one should not expect an oscillatory behavior of Hq . However, according to 
our calculations of the correlators with the MC generators JETSET 7.3 ||21[], ARIADNE 
4.4 [|l(e+e-, =91 GeV) and PYTHIA 5.5 EKe-p, =314 GeV) Hq has, nev- 
ertheless, an oscillating form as shown in Fig. 4. An explanation of such a phenomenon 
can be found immediately if one recall two facts: 1) each MC generator takes a special 
care about both the local (in the course of parton splitting) and global energy-momentum 
conservation in the collision; 2) the finite energy of collisions is the physical origin of large 
0{as) corrections Thus, the LLA in conjunction with the energy- momentum conser- 
vation in the MC models imitate in some part the higher order corrections leading to the 
oscillation of Hq . The question arises, though, how much of the higher order corrections 
are accounted for? 
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Fig. 5. The ratio Hq calculated for NBD with (n)=9.22, k=17.2A. The solid line are 
for NBD truncated at ntr and the dished lines are for NBD without truncation. 



5 Phenomenological examples 



The conclusions from the previous section can be confirmed by the following arguments 
||22|| . Formally, according to (7) NBD has an infinite "tail" at finite collision energy (finite 
(n)). This results in positive Kg and monotone declining Hq (8). On the other hand, 
an infinite "tail" of MD is possible only for production of massless particles or through 
neglecting energy conservation during the reaction. Taking into account these factors 
leads to a truncation of the MD "tail" at some finite multiplicity ntr{s). As a result, Hq 
calculated for the truncated NBD oscillates around the curve with alternating sign. 
The amplitude of the oscillation tends to zero quickly 

^(NBD) ^pig_ 5) 



oo and if^**") tends to 



The same behaviour of Hq has been found for the truncated Poisson 
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distribution (PD). 

Another example of the behaviour of Hq is in soft pp colhsions at SppS and Tevatron 
energies calculated in ||2^ in the framework of the Dual Parton Model [24|. It was found 
23[| the properties of Hg ( amplitude of the oscillation, positions of minima and maxima) 



are very sensitive to the number of cutted Pomerons accounted for in the calculation. 



6 What can be done at HERA ? 



In high energy reactions used in the study of the oscillations of Hq |[T^ only DIS data are 



missing. New data from the ep collider HERA will be able to rectify this situation. The 
invariant mass W of the hadronic final state in DIS at HERA extends, with significant 
cross sections, to the phase space limit {y/s =314 GeV). This circumstance allows us to 
formulate several problems related to properties of MD which can be studied with the HI 
and ZEUS detectors: 

1. Detailed study of MD as a function of z = n/ (n) over the whole kinematical region 
of W. Does the KNO scaling violated at large W ? 

2. High precision measurement of the ration Hq = Kq/Fq of cumulant and factorial 
moments both for the full phase space and restricted rapidity windows, for events with 
1+1, 1+2, ... jets, etc. Does Hq as a function of the order q shows an oscillations around 
Hq = 07 If so, confronting the data with predictions of the MC models we would learn 
more about the higher order effects implemented in these MC models. 

3. Measurements of Hq at different W will shed more light on the problem how the 
finite energy effects influence the Hq shape. 

To conclude, energy-momentum conservation plays a very important role in the correct 
description of the multiplicity distributions , Fq , Kq and Hq in the framework of QCD 
and different phenomenological models. The ratio Hq is extremely sensitive to the length 
of the MD "tail". In perturbative QCD the behaviour of the MD "tail" is controlled by 
higher order corrections while for phenomenological approaches (NBD, PD etc.) the finite 
energy effects have to be accounted for by truncating the MD "tail". 
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